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Abstract 

Among all neutrino mixing parameters, the atmospheric neutrino mixing angle #23 
introduces the strongest variation on the flux ratios of ultra high energy neutrinos. 
We investigate the potential of these flux ratio measurements at neutrino telescopes 
to constrain #23- We consider astrophysical neutrinos originating from pion, muon- 
damped and neutron sources and make a comparative study of their sensitivity reach 
to 023- It is found that neutron sources are most favorable for testing deviations from 
maximal #23- Using a x 2 analysis, we show in particular the power of combining (i) 
different flux ratios from the same type of source, and also (ii) combining flux ratios 
from different astrophysical sources. We include in our analysis "impure" sources, 
i.e., deviations from the usually assumed initial (1 : 2 : 0), (0 : 1 : 0) or (1 : : 0) flux 
compositions. 
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1 Introduction 



Dedicated facilities spanning km 2 of area for detecting ultra high energy neutrino com- 
ing from astrophysical sources are under construction (IceCube [1,2]) or consideration 
(KM3Net [3]). Motivated by this, a significant number of papers has been devoted in 
recent years to the phenomenon of neutrino mixing of high-energy astrophysical neutri- 
nos [4-22]. That neutrinos are massive and therefore mix has been proved beyond any 
doubt by observations of neutrinos coming from the Sun [23], atmosphere [24], reactors [25] 
and accelerators [26,27]. The most recent limits on the mass-squared differences and the 
mixing angles can be found in [28]. The mass splitting associated with solar neutrino 
oscillations is Am^ ~ 7.6 x 10~ 5 eV 2 , while that associated with atmospheric neutri- 
nos is jAm^l ~ 2.5 x 10~ 3 eV 2 . Because the oscillation lengths corresponding to these 
mass-squared differences are much smaller than astrophysical distances, the oscillations 
get averaged out for the astrophysical neutrinos. However, the non-trivial flavor mixing 
in the lepton sector still modifies the neutrinos fluxes on their way from the source to the 
detector. This opens up the possibility to obtain information on the mixing angles and the 
CP phase of the PMNS matrix. This information would be complementary to the already 
impressive existing and expected data from current and future experiments devoted purely 
to neutrino oscillations. 

Ultra high energy (UHE) neutrinos are expected to come from decay of pions, muons 
and/or neutrons. Therefore, even though the absolute numbers of UHE neutrinos are 
uncertain by a huge amount, the relative proportions of the initial flux compositions $g : 
<£>° : <£>° are known. Here $° with a — e, n, r is the initial flux of a neutrino with flavor 
a. This ratio is (1 : 2 : 0), (0 : 1 : 0) and (1:0:0) for pion, muon-damped and neutron 
sources, respectively. Hence, working with the flux ratios of different flavors is considered 
to be much less model dependent than working with absolute fluxes. One complication 
which still arises is that in general one expects corrections [8,19,30] to the usually assumed 
"pure" initial flux compositions. For instance, instead of ($g : $° : $°) = (1:2:0), one 
might initially have (1 : 2 (1 — () : 0), with ( being around = 0.1 [30]. Not taking such 
deviations into account can lead to wrong conclusions about the neutrino parameters [22] . 

Currently the mixing angle 6*23 is one of the least known besides #13 , the CP phase and 
the neutrino mass hierarchy. In fact, the quantity sin 2 #23 is uncertain by about ±33% 
at 3er. Information on this mixing angle and its deviation from maximality can be tested 
in future atmospheric neutrino experiments [31,32]. Using the zenith angle dependence 
of the muon events in atmospheric neutrino data, one would be able to restrict sin 2 #23 
to within ±24% with 1.84 MTy (megaton-year) data in water Cerenkov detectors and to 
within ±30% with 250 kTy data in large magnetized iron calorimeters [33]. Atmospheric 
neutrino data could also be used very effectively to give us the "octant" (whether #23 < vr/4 
or > 7r/4) of #23, if indeed #23 turns out to be non-maximal. The sub-GeV electron 
events in water Cerenkov detectors carry information on the #23 octant through the Am^- 
driven sub-dominant oscillations [31]. The multi-GeV muon events in large magnetized 
iron calorimeters have sensitivity to Earth matter effects which in turn depend on 623 and 
its octant [32]. Long baseline experiments would also give very good sensitivity to sin 2 #23: 
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the combined 5 year data from MINOS, ICARUS, OPERA, T2K and NOz/A is expected 
to constrain sin 2 #23 to within ±20% around its maximal value [34]. 

In this letter we will focus on the possibility to constrain the atmospheric neutrino 
mixing angle #23 from measurements of flux ratios at neutrino telescopes. After arguing 
that their dependence on this angle is strongest among the neutrino mixing parameters, we 
discuss its possible constraints using a x 2 analysis. We quantify that combining different 
sources leads to stronger constraints on #23 • If different flux ratios from different neutrino 
sources are combined (those are pion, muon-damped and neutron sources), even better 
constraints are possible. Neutron beam sources turn out to be the most interesting ones. 
We include the possibility of "impure sources" in our analysis and investigate their impact 
for the first time in a statistical analysis. 

We begin by discussing the flux ratios of UHE at the neutrino telescopes in Section 
2. In Section 3 we introduce our \ 2 function and use it to give prospective bounds on #23 
using the flux ratios. We end in Section 4 with a summary of our results and conclusions. 



2 Neutrino Mixing and Neutrino Telescopes 

Astrophysical sources will generate fluxes of electron, muon and tau neutrinos, denoted by 
<£>e, and respectively. As a consequence of non-trivial neutrino mixing, it is not this 
initial flux composition which arrives at terrestrial detectors. In fact, what is measurable 
is given by 
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and U is the lepton mixing matrix. The flavor mixing matrix P comprises of the individual 
P a p elements. The current best-fit values as well as the allowed la and 3a ranges of the 
oscillation parameters are [28] : 
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These mixing angles can be related to elements of the PMNS mixing matrix via 
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where 



cos%, Sij = sin 6ij. The CP phase S is unknown. Because the mass-squared 



differences drop out of the mixing probabilities, and solar neutrino mixing is neither max- 
imal, zero nor tt/2, no transition probability P a p with a 7^ (3 is zero and no survival 
probability P aa is one. Consequently, high energy astrophysical neutrinos will always mix. 
To be precise, at la and 3a the entries of the flavor conversion matrix P a p are 



P 
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0.34-^0.44 0.35 
0.35 
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As is obvious from Eq. 
bimaximal mixing [35], 



(at la) , 
(at 3a) . 

3]), a good zeroth order description of neutrino mixing is tri- 
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Therefore, it proves in particular useful to expand in terms o£] 
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For simplicity, we fix from now on for analytical considerations sin 2 # 12 to |. This is a valid 
approximation for obtaining simple analytical expressions, but we stress that the numerical 
results to be presented in this paper are obtained with the exact expressions. The result 
for the flavor mixing matrix is 
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where the universal first [12,15] and second [21,22] order correction terms are 

1 
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A 2 = -(2cos 2 5|t/ e3 | 2 + 7e 2 - V2 cos5\U e3 \e 
9 ^ 



(9) 



We have omitted here small and usually negligible terms in the expansion which only 
contain |?7 e 3| 2 , see Refs. [21,22]. Note that A can be written as | (2y / 2 cos S \ U e s\ — e) 2 +3 e 2 



1 An expansion up to second order around zero \U e s\, maximal #23 and sin 2 #12 = 5 can be found in 
Refs. [21,22]. 
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and therefore is positive semi-definite. In addition, A turns out to be often larger than 
the first order term A [21,22]. To be quantitative, 

at la : -0.043 < A < 0.069 , A 2 < 0.061 , 

" " _ 2 ~ (10) 

at 3(7 : -0.104 < A < 0.117 , A < 0.179 . 

The dependence of A and A on 812 is very weak [21,22]. From the expressions for A 

— 2 

and A it is clear that their dependence on the atmospheric neutrino mixing angle 823 is 
stronger than on cos 5 \ U e s\. It is the large prefactor in front of e and e 2 which is the reason 
for this behavior. In addition, these prefactors are larger than the ones related to \U e s\. 
Hence, the dependence on |{7 e3 | is weaker and smeared by the additional dependence on 
cos 5. This is why we are interested here in effects of deviations from maximal atmospheric 
neutrino mixing. 

The neutrino sources we assume are pion, muon-damped [36] and neutron [37] sources, 
with initial flux compositions of 

C (1:2(1-C):0) pion, 
($0 : $0 : $°) = ) (77 : 1 : 0) muon-damped , (11) 
[ (1 : 77 : 0) neutron . 

By introducing small (,?j<1 [22], we have allowed here for impurities in the initial com- 
positions, which can be expected on general grounds [8,19,30]. 

Turning to the observable flavor ratios [38,39], the most frequently considered is the 
ratio of muon neutrinos to all other flavors: 

T= ^ = -^. (12) 

Using again sin 2 6*12 = | and expanding in terms of the small parameters | U e 3 \ , e — | — 823 
as well as ( or rj one finds [22] 

|(1-A + A 2 -|C) , pion (1:2(1-C):0) , 
T ~ \ X_A + iA 2 -|r], muon-damped (77 : 1 : 0) , (13) 
I + A + I r\ , neutron beam (1 : rj : 0) . 

— 2 

The second order correction A appears only in P w , P^ T and P TT , and therefore does 
not affect T for neutron sources. Hence, the dependence on sin 2 6*23 can be described to 
an excellent approximation as quadratic for pion and muon-damped sources, but linear for 
neutron sources. We show in the left panel of Fig. [T]the ratio T as a function of sin 2 823 and 
sin^ for the three pure sources (i.e., ( = rj = 0). Plots are shown for fixed sin 2 6 , 12 = 0.32 
and three values of S. The stronger dependence on sin 2 #23 is clearly seen. Fig. [2]focusses on 
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the dependence of the flux ratios on sin 2 623- We have taken the range sin 2 #12 = 0.32±0.02, 

sin 2 # 13 < 0.005 and < 5 < 2n. We have also marked the current la and 3a allowed 

ranges of sin 2 #23 for comparison. We reiterate that these plots are for ( = r\ = 0. Note 

that because of larger prefactors in front of ( and r\ in Eq. (flBl) . impurities are expected to 

have stronger impact for muon-damped and neutron sources. Finally, we stress that the 

2 

correction factors A and A , which can both be (9(0.1), are added to a number around 1 
for pion sources, but to a number around 0.4 for muon-damped and 0.2 for neutron sources. 
Hence, neutron sources are expected to be the most useful to constrain #23 • 
Another ratio which can be considered is the ratio of electron to tau neutrinos, 



For the three neutrino sources one finds [22] 



Comparing these expressions with the ones for T in Eq. fflBl . we note that in case of pion 

— 2 

and muon-damped sources the first and second order correction terms A and A are added 

in R, whereas they are subtracted in T. Recalling that e and e 2 appear with equal sign 
2 

in A and A , we expect stronger dependence on (9 23 in R than in T. This is illustrated in 
the right panels of Figs. [Hand El where we display R as a function of sin^ and sin 2 #23 . 
Again, the dependence on sin 2 623 is basically quadratic for pion and muon-damped sources, 
whereas it is linear for neutron beams. The effect of impurity will be in general larger for 
R, because of the larger prefactors in front of ( and rj. Note also that the parameters ( 
and 7/ appear in R with opposite sign compared to T. 

3 Statistical Analysis 

To statistically investigate the prospects of constraining 6*23 with neutrino telescopes, we 
turn to a x 2 analysis. First we discuss the case in which only T is measured. Let us define 
the x 2 function to be minimized as: 



where T t h and T exp are the theoretically predicted and experimentally measured T ratios, 
respectively, and <7t cxp is the la uncertainty of the measured value T exp in the neutrino 
telescope. 




pion (1:2(1-C):0), 
muon-damped (rj : 1 : 0) 
neutron beam (1 : 77 : 0) 



(15) 




(16) 
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We take into account the possible uncertainty of the initial flux composition by con- 
structing different x 2 functions for different initial flavor compositions and then define the 
allowed region of sin 2 823 as the maximum range obtained by combining the different func- 
tions. The errors on T for the different sources are expected to lie in the range 15-25% 
after a decade of of running of the experiment. It is expected that the error on pion sources 
will be less than those on muon-damped and neutron sources. Therefore, we present most 
of our results for fixed assumed errors on the flux ratio T as 10% for pion, 20% for muon- 
damped and 15% for neutron sources. In what regards R, the assumed errors are 15% for 
pion, 25% for muon-damped and 20% for neutron sources. In general one expects larger 
errors for R than on T: muon neutrinos can be identified at neutrino telescopes for energies 
greater than 100 GeV through characteristic muon tracks. Electron and tau neutrinos, can 
be identified through showers and effects like "double bang" events, respectively, only for 
higher energies, above 10 6 GeV. Recall that muon-damping is expected to happen in a 
generic pion source for a particular and limited energy range, in which the muon is ab- 
sorbed before decaying [36]. Hence, less statistics for muon-damped sources with respect 
to pion sources is expected. Concerning neutron sources, we have chosen an error greater 
than the one for pion sources, because they are expected to be characterized in an energy 
range covered by atmospheric neutrino background [37]. Therefore the systematic error 
in this case might be greater than the one for pion sources. Hence, we have chosen a 
hierarchy among the errors for the different types of sources. However, the values for the 
errors taken above are just one choice for the errors, which could be different from what 
has been chosen above. To take into this fact and to show the impact of the errors on our 
analysis, we will also show the results as a function of the error on T. 

We generate the prospective data T exp for two sets of mixing angles which we assume 
as "true" and minimize the \ 2 function to obtain bounds on the measured 823 at neutrino 
telescopes. In the fit we allow 6 12 , 8 13 and 5 to take any value in their physically allowed 
ranges. The second term in Eq. fflBI) takes into account the "priors" on the mixing angles 
9i2 and #13, on which we expect better constraints by the time we get the data on UHE 
neutrinos at neutrino telescopes. For the la uncertainty on 6* 12 we use the range given in 
Eq. ([3D, while for # 13 we use the following upper bound: 

sin 2 20 13 < 0.03 at 90% C.L. , (17) 

which corresponds to the absence of a signal in the Double CHOOZ experiment after three 
years of operation with both detectors [40]. We have considered in our analysis this specific 
case and we have studied the consequences that can be inferred on the atmospheric mixing 
angle. This situation has be considered as an optimistic scenario. Indeed, larger values of 
813 than the ones we have chosen would worsen the sensitivity on 823- 

We generate T exp and show results for two sets of true values for the mixing parameters: 

1. sin 2 ^12 = 1/3 , sin 2 ^13 = and sin 2 #23 = 1/2 (Scenario TBM) , 

2. sin 2 12 = 0.32 , sin 2 8 l3 = and sin 2 823 = 0.6 (Scenario 2) . 
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Scenario TBM gives T exp (and R exp ) corresponding to their "tri-bimaximal values", i.e., the 
zeroth order terms from Eqs. ({TBI and (|T5|) . Scenario 2 gives T cxp = 0.36, 0.45, and 0.17 
respectively for the pion, muon-damped and neutron sources. The corresponding values of 
-Rexp are given as 0.90, 0.45, and 2.16 respectively. Scenario TBM corresponds to the case 
where the true value of 623 is maximal, while Scenario 2 exemplifies a situation where the 
true value of 823 is non-maximal and > 7r/4. We checked that the results for 823 < 7r/4 
share the same features as the ones for 823 > 7r/4. 

We show the results of our \ 2 fit using T only in Fig. [33 The left panels show results 
for Scenario TBM while the right panels are for Scenario 2. The upper, middle and lower 
panels are for pion, muon-damped and neutron sources, respectively. In each panel the 
current lcr and 3a limits on sin 2 823 are given. We have allowed for impure sources by 
choosing initial flux compositions of (1 : 2 : 0), (0 : 1 : 0) and (1:0: 0), (1 : 1.9 : 0), 
(0.05 : 1 : 0) and (1 : 0.05 : 0), as well as (1 : 1.8 : 0), (0.1 : 1 : 0) and (1 : 0.1 : 0) 
for pion, muon-damped and neutron sources, respectively. Interestingly, impurity affects 
the neutron sources more strongly, and the impact is slightly larger on the "dark-side" of 
sin 2 #23, i.e., for 823 > tt/4. The "mexican hat" shape for pion and muon-damped sources 
is easy to understand by looking at the sin 2 823 dependence of T in Fig. [2j Therefore, for 
pion and muon-damped sources we get, in general, a two-fold degeneracy for every value 
of T. For Scenario 2 we get slightly larger values for T and are hence no longer at the 
bottom of the parabolic shape. Therefore, the x 2 minima lie at lower and larger values of 
sin 2 823 and the two possible values of sin 2 823 become further separated. 

From Fig. [3] we can infer that considering only the pion or muon-damped flavor ratio 
will bring limited information on sin 2 823 . Indeed for these two sources and for Scenario 
TBM, we get an allowed lcr range that is much bigger than the current la range. With 
neutron sources, we get an allowed la range for sin 2 823 from 0.42 to 0.64, which is a bit 
smaller than the current 3a range. For Scenario 2 (true value sin 2 823 > 1/2) we would 
not be able to corroborate it with either the pion or muon-damped sources, in particular 
because we obtain two allowed regions for sin 2 823- With neutron sources, and their linear 
dependence on sin 2 823 , we expect to get only one allowed zone. The plot confirms this 
and in addition we see that in this case maximal mixing is ruled out by neutrino telescope 
data alone. To be more precise, from measuring T with neutron sources we could constrain 
sin 2 823 to be greater than 0.54 at the la level and greater than 0.5 at 90% C.L. 

Since these conclusions depend on the errors considered for the flavor ratios, we have 
generalized our analysis and studied the allowed range of sin 823 as a function of the error. 
In Fig. H] we display the results. The upper 4 panels are for data generated for the TBM 
case while the lower 2 are for the neutron sources only and data corresponding to Scenario 
2. For TBM in data and with pion and muon-damped sources, we can infer that even for 
an extremely small error of the order of 5%, the allowed range at 90% C.L. is bigger than 
the current ones. Neutron sources can impose an upper limit only for very small errors. A 
lower limit can be imposed if the error on the flavor ratio is smaller than 18%. If Scenario 
2 was true, then neutron sources could rule out maximal mixing and give information on 
the octant of 823, and Fig. [5] shows that the T ratio is better suited than R for this purpose. 

To improve the sensitivity of neutrino telescopes to 823, one should combine measure- 
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ments of T from different sources. In Fig. [5] we have considered the combination of pion 
and muon-damped sources (upper panels), of pion and neutron sources (middle panels) and 
of all the three different type of sources (lower panels). The \ 2 function is, for example, 
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(18) 



when measurements from pion and muon-damped sources are combined. The left panels 
show results for the Scenario TBM, while the right panels showcase Scenario 2. It should 
be noted that the two-fold degeneracy with respect to sin 2 9 2 3 is hardly present anymore 
when two or more T ratios are combined. Three different cases are plotted in the figures. 
With "pure" we refer to initial flavor compositions (1 : 2 : 0), (0 : 1 : 0) and (1:0: 0), 
with "impure A" we refer to (1 : 1.9 : 0), (0.05 : 1 : 0) and (1 : 0.05 : 0), while for 
"impure B" we refer to (1 : 1.8 : 0), (0.1 : 1 : 0) and (1 : 0.1 : 0). Again, the effect 
of impurity is stronger for sin 2 # 23 > 1/2. In general, impurity shifts the \ 2 minimum to 
larger values of sin 2 #23 • Analyzing the plots, we can conclude that the combinations of 
pion and muon-damped sources will not give us precise information on sin 2 #23 and that 
only the combination with neutron sources improves substantially the sensitivity. Indeed, 
we can notice that the combined constraint from pion and neutron sources looks basically 
identical to the result of neutron sources alone (see Fig. [3D- Moreover, if all three sources 
are combined, i.e., if muon-damped sources are added, hardly any improvement is seen. 
For Scenario 2, we could determine the octant of #23 at the la level in case we combine T n 
and T n or all the three flavor ratios. We can, therefore, conclude that the importance of 
the neutron source is crucial to provide a good sensitivity to sin 2 # 23 . 

Better constraints on #23 are obtained when the ratios T and R are combined. The % 2 
function is now the same as in Eq. (1161) with an appropriate term including R added. The 
result of the minimization is shown in Fig. [61 and the improvement with respect to the 
other plots is obvious. The effect of impurity is however stronger. This is understandable 
from the approximate expressions of T and R in Eqs. fTl3l[T5|) . in which the impurity factors 
( and 77 have more sizable prefactors in R. 

To summarize, once one combines T and R, the x 2 curves become less broad, but the 
shift due to impurity becomes more sizable. The best constraint is obtained by combining 
T and R for neutron sources. Indeed, Fig. [6] shows that the la range for sin 2 #23 is of the 
same order as the current one. For Scenario 2 we can exclude maximal mixing at roughly 
90% C.L. level. This conclusion depends also on the error considered for R. Varying it we 
show in Fig. H] the allowed range of sin 2 # 23 if only R was measured. 

The most optimistic scenario would occur when both ratios can be measured for all 
three sources. We checked that the outcome is basically identical to combining T and R for 
neutron sources, showing once again that these sources are best suited for 9 2 3 constraints. 

We stress that throughout this study we have tried to extract information on only the 
mixing angle (9 23 and assumed that the kind of astrophysical source is known. Note that 
the values for the T and R ratios of different sources are well separated even if the errors 
on the flux ratios are bigger than the ones we are considering. In addition, the neutron 
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beams, for which the sensitivity to sin 2 #23 is stronger, have the advantage that the ratio T 
is rather small, and R is very large. This will allow more easily to distinguish one source 
from the others. A detailed analysis on the problematics related to the determination of 
the type of source is beyond the task of this work and will be published separately. 

4 Summary and Conclusions 

Neutrino mixing affects the flux ratios of ultra high energy neutrinos arriving on earth. 
The exact values of the ratios are determined by the values of the mixing parameters. A 
series of papers have looked into the impact of standard and non-standard properties of 
neutrinos on the flux ratios. In this paper we argued that within the framework of the 
standard picture with three stable neutrinos, the mixing angle #23 has the maximum effect 
on the flux ratios. We performed a statistical test to ascertain quantitatively the extent 
to which this mixing angle can be constrained by measuring the flux ratios of ultra high 
energy neutrinos. 

We defined two kind of flux ratios, T = $ M /$ tot and R = $ e /$ r , where <3> tot = $ e + 
®n + $r- We showed the dependence of these ratios on the mixing angles 6> 13 and #23 for 
three kinds of sources of ultra high energy neutrinos - pion, muon-damped and neutron 
sources. Their dependence on #23 is largest, and from analytical considerations we gave an 
idea of the extent to which the flux ratios could be used to constrain it. 

We defined a \ 2 function to quantify the extent to which #23 can be constrained by 
neutrino telescopes. We assumed a simplistic approach where we worked with the flux 
ratios themselves but we are aware that in a realistic analysis one should work with ratios 
of actual number of events, taking into account the detector response and efficiencies. 
However, we have chosen to work in a simplified set-up because the purpose of this paper 
was to make a ballpark estimate of the sensitivity of neutrino telescopes to #23- For the first 
time we included "impure" initial neutrino fluxes in a statistical analysis by introducing 
two variables, ( and rj, which parameterize deviations from the initial flux compositions 
(1 : 2 : 0), (0 : 1 : 0) or (1 : : 0) for the pion, muon-damped and neutron sources. 

We performed the statistical analysis using T from one given source at a time and 
conclude that the best #23 sensitivity comes from neutron sources. We showed that an 
error of less than 20 % is necessary to obtain results comparable to oscillation experiments. 
We presented results by combining pion and neutron sources, muon-damped and neutron 
sources, and finally all three taken together. The bound improves mainly because neutron 
sources have a much better handle on #23 • We next combined T and R measurements at the 
neutrino telescopes. Adding the information of R improves the #23 bound substantially, 
but increases somewhat the impact of impure sources. In particular, we note that the 
combination of T and R for neutron sources could give us bounds on #23 which are much 
better than the one we have currently. In fact, measuring T and R for neutron sources could 
be comparable to the bounds on #23 we expect from future long baseline and atmospheric 
neutrino experiments, depending of course on the uncertainty on the measurement of the 
flux ratios at the neutrino telescopes. We performed the statistical test on #23 for true 
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sin 2 #23 = 0.5 and sin 2 #23 = 0.6. For the latter case we checked if it would be possible to 
establish the right octant of 6*23. 

In conclusion, in favorable but not unrealistic situations we can indeed expect useful and 
complementary limits on #23 which are comparable to the ones from dedicated oscillation 
experiments. 
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Figure 1: The ratios T = $ M /$tot (left panel) and R = $ e /$ r (right panel) for the three 
sources under discussion as a function of sin#i 3 and sin 2 6*23. The mixing angle 9 12 has been 
fixed to its best-fit value, i.e., sin 2 6*12 = 0.32. 
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Figure 2: The ratios T = $ M /$ t ot (left panel) and R = $ e /$ r (right panel) for the three 
sources as a function of sin 2 # 2 3 when the ranges of the mixing parameters are sin 2 # 12 = 
0.32 ±0.02 and sin 2 # 13 < 0.005. The current la and 3a ranges of sin 2 623 are also indicated. 
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Figure 3: Constraints on sin 2 623 from a x 2 analysis obtained by measuring the ratio 
T = $ M /$ to t for the three different neutrino sources. The left panel assumes the "tri- 
bimaximal values" T — |, ^ and | as experimental values, while the right panel is for 
"true values" of the mixing angles given by Scenario 2. The current la and 3a ranges of 
sin 2 #23 are also indicated. 
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Figure 4: Allowed range of sin $23 3s a function of the error on the flavor ratio. The dark 
areas represent the la range, while the lighter areas are at 90% CL. The upper 4 plots 
are for scenario TBM, the 2 lower ones for scenario 2, where we only plot neutron sources. 
The current la and 3<r ranges of sin 2 $23 are also indicated. 
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Figure 5: Same as Fig. [31 but for combined measurements of T from different sources. 
With "pure" we refer to the initial flavor composition (1 : 2 : 0), (0 : 1 : 0) and (1:0:0), 
with "impure A" we refer to (1 : 1.9 : 0), (0.05 : 1 : 0) and (1 : 0.05 : 0), while for "impure 
B" we refer to (1 : 1.8 : 0), (0.1 : 1 : 0) and (1 : 0.1 : 0). 



